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Abstract 



We consider QCD radiative corrections to WW pair production in hadron collisions. 
We perform a calculation that consistently combines next-to-leading order predictions 
with soft-gluon resummation valid at small transverse momenta p^ w of the WW pair. 
We present results for the p^ w distribution at the LHC up to (almost) next-to-next- 
to-leading logarithmic accuracy, and study the effect of resummation on the charged- 
lepton distributions. Soft-gluon effects are typically mild, but they can be strongly 
enhanced when hard cuts are applied. The relevant distributions are generally well 
described by the MC@NLO event generator. 
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1 Introduction 



Vector boson pair production in hadron collisions is important for physics within and beyond the 
Standard Model (SM). First of all, this process can be used to measure the vector boson trilinear 
couplings. Any deviation from the pattern predicted by SU(2) (g> U(l) gauge invariance would 
be a signal of new physics. The Tevatron collaborations are currently measuring WW, WZ and 
W'y cross sections at invariant masses larger than those probed at LEP2, setting limits on the 
corresponding anomalous couplings [lj. 

Furthermore, vector boson pairs are an important background for new physics searches. If a 
heavy enough Higgs boson does exist, it will decay with large branching ratios in WW and ZZ 
pairs, whereas charged Higgs bosons from non standard Higgs sectors may decay into WZ final 
states [2]. Typical signals of supersymmetry, e.g. three charged leptons plus missing energy, find 
an important background in WZ and Wj production [3]. 

In this context, WW production has received considerable attention. This process provides 
an irreducible background to the Higgs boson search in the H — > WW —>■ lulu channel, which is 
the most important when the mass Mh of the Higgs boson is in the range 155 < Mh ^170 GeV. 
Since two neutrinos are present in the final state, no invariant-mass peak can be reconstructed. 
Fortunately, the angular correlations among the charged leptons suggest a good discrimination of 
the signal over the background [I] , provided the background can be reliably extrapolated into the 
signal region. It is thus essential to have good control on the SM prediction for the WW cross 
section as well as for the associated distributions. 

QCD corrections to WW production at next-to-leading-order (NLO) have been computed 
more than 10 years ago [5JE] and enhance the cross section by about 40% at LHC energies. These 
calculations were done with the traditional method of evaluating directly the relevant squared 
amplitudes. As a consequence, the W's polarization was summed over and spin correlations were 
not taken into account. 

More recent NLO calculations exist that, using the one-loop helicity amplitudes computed in 
Ref. [7], fully take into account spin correlations [8]. The general purpose NLO program MCFM 
[H] includes in addition single- resonant contributions neglected in the calculation of Ref. [SJ . Very 
recently, the effect of the potentially large next-to-next-to-leading (NNLO) contribution from the 
one-loop gg — > WW — > lulu diagrams has also been studied [TUJ [TT] . 

The fixed-order NLO calculations provide a reliable estimate of WW cross sections and distri- 
butions as long as the scales involved in the process are all of the same order. When the transverse 
momentum of the WW pair p^f w is much smaller than its invariant mass M^yiy the validity of the 
fixed-order expansion may be spoiled since the coefficients of the perturbative expansion can be 
enhanced by powers of the large logarithmic terms, ln n My/w / Pt W ■ This is certainly the case for 
the pj? w spectrum, which, when evaluated at fixed order, is even divergent as pj? w — > 0, and thus 
requires an all-order resummation of the logarithmically enhanced terms. Resummation effects, 
however, can be visible also in other observables, making it important to study them in detail. 

The way to perform transverse- momentum resummation is known [12]- [2D]: to correctly im- 
plement momentum conservation in the transverse plane the resummation has to be performed in 
impact parameter 6-space, b being the variable conjugated to pt through a Fourier transformation. 
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In the case of vector boson pair production, transverse-momentum resummation has been applied 
to 77 [21] and ZZ [22] pair production. To our knowledge, the case of WW pair production has 
not been considered yet. 

In the present paper we report on an implementation of 6-space resummation to WW pro- 
duction in hadron collisions. We use the helicity amplitudes of Ref. [7] and work in the narrow 
width approximation (i.e. we only consider double-resonant contributions), but fully include the 
decay of the W*s, keeping track of their polarization in the leptonic decay. In the large p^ w 
region we use LO perturbation theory parton); in the region p^ w <C Mww the large 

logarithmic contributions are resummed to (almost) next-to-next-to-leading logarithmic (NNLL) 
[2"3| 121] accuracy. Our results have thus uniform NLO accuracy but consistently include the effect 
of resummation in the region p^ w <C M ww . 

To perform the resummation we use the formalism of Refs. [251 [26] . In this approach, the 
resummation is achieved at the level of the partonic cross section and the large logarithmic contri- 
butions are exponentiated in a process-independent manner, being constrained to give vanishing 
contribution to the total cross section. 

Besides computing transverse momentum spectra of the WW pair, we study the effect of 
resummation on the leptonic distributions, when different cuts are applied. We also compare our 
results with those obtained at NLO and with the ones from the general purpose event generator 
MC@NLO [27], which, in its latest release [28j, partially includes the effect of spin correlations in 
the Ws decay. 

The paper is organized as follows. In Sect. [2] we discuss the application of 6-space resummation 
to WW production. In Sect. [3] we present our numerical results and in Sect. 0] we draw our 
conclusions. 

Preliminary results of this work for the p^ w spectrum were used in the study of Ref. [29J . 



2 Transverse- momentum resummation for WW pair pro- 
duction 

In this Section we apply 6-space resummation to the production of WW pairs in hadron collisions. 
The resummation formalism we use is completely general, and, as discussed in detail in Ref. [26], 
can be applied to a generic process in which a system of non strongly interacting particles of high 
mass M is produced in hadronic collisions. In the case of WW production, the mass M is the 
invariant mass of the WW system. 

The resummation is performed at the level of the partonic cross section, which is decomposed 

as: 

,~ 7 * (res.) r~ (fin.) 

aVWWab _ dV W Wab , d(J WWab /,x 

dM 2 dp 2 T dM 2 dp 2 T dM 2 dp 2 T ' 1 ' 

The first term on the right hand side, da^^ ab , contains all the logarithmically enhanced con- 
tributions at small pj>, and has to be evaluated by resumming them to all orders in a s . The 
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(fin.) 
WWc 

perturbation theory. 



second term, da WV y ab , is free of such contributions, and can thus be evaluated at fixed order in 



The resummed component da^ w \ b can be expressed as 

^^(PT,M,s;a s (^ R )^ 2 R ^l) = — J db - J (bp T ) W% w (b, M, s; e*s(/i|), /4 /4) , (2) 

where Jq(x) is the 0-order Bessel function, fi R (hf) is the renormalization (factorization) scale 
and s is the partonic centre-of-mass energy. By taking the iV-moments of W with respect to the 
variable z = M 2 /s at fixed M the resummation structure of W^ b N can indeed be organized in 
exponential formQ 

W% w (b, M; a s (/4), /4 /4) = H% w (M, a s (/4); M 2 //4 M 2 //x|, M 2 /Q 2 ) 

x exp{g N (a s (fi 2 R ), L; M 2 //4 M 2 /Q 2 )} , (3) 

were we have defined the logarithmic expansion parameter L as 

2 b 2 

L^y~ (4) 

and the coefficient b = 2e _7B (7^ = 0.5772... is the Euler number) has a kinematical origin. 
The scale Q appearing in Eqs. ([2J HJ), named resummation scale in Ref. parametrizes the 
arbitrariness in the resummation procedure, and has to be chosen of the order of the hard scale 
M. Variations of Q around M can give an idea of the size of yet uncalculated higher-order 
logarithmic contributions. The function 7ijy w does not depend on the impact parameter b and it 
includes all the perturbative terms that behave as constants as b — > 00. It can thus be expanded 
in powers of «s = «S : 

H% W (M, as; M 2 //4, M 2 //4, M 2 /Q 2 ) = a$ w (<*s, M) [l + f < W (1) (M 2 /^, M 2 /,a|, M 2 /Q 2 ) 



+ p) < W(2) (M 2 /,a 2 „M 2 / / a|,M 2 /Q 2 ) + ...] (5) 



where u w ^ w is the Born partonic cross section. The exponent Qn includes the complete dependence 
on b and, in particular, it contains all the terms that order-by-order in as are logarithmically 
divergent as b — > 00. The logarithmic expansion of Qn reads 

g N (a s (fi 2 R ),L;M 2 /ii 2 R ,M 2 /Q 2 ) = Lg® (a s L) + g%\a s L; M 2 / /4 M 2 /Q 2 ) 

+ %£W, M 2 //4 M 2 /Q 2 ) + . . . (6) 

7T 

where the term L collects the LL contributions, the function includes the NLL Contribu- 
te) 

tions, g N controls the NNLL terms and so forth. 

In the implementation of Eq. ([3]) the resummation of the large logarithmic contributions affects 
not only the small-pr region {jpT <C M), but also the region of large pt {pt ~ M). This can 



^Here, to simplify the notation, flavour indices are understood, or in other words, we limit ourselves to discussing 
the flavour non-singlet contribution. A complete discussion of the exponentiation structure in the general case can 
be found in Appendix A of Ref. . 
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be easily understood by observing that the logarithmic expansion parameter L is divergent as 
b —>■ 0. To reduce the impact of unjustified higher-order contributions in the large-pr region, the 
logarithmic variable L in Eq. (j3J) is replaced as 



M^ + l). (7) 



q% 2 

b 2 

The variables L and L are equivalent when Qb 3> 1, but they lead to a different behaviour of 
the form factor at small values of b (i.e. large values of pr)- When Qb <C 1 in fact, L — > and 
expj^iv} — > 1. The replacement in Eq. ([7]) has thus a twofold consequence: it reduces the impact 
of resummation at large values of pt, and it allows us to recover the corresponding fixed-order 
cross section upon integration over px- 

Another important property of the formalism of Ref. [26] is that the process dependence (as 
well as the factorization scale and scheme dependence) is fully encoded in the hard function 7i ww . 
In other words, the functions are universal: they depend only on the channel in which the 
process occurs at Born level, (qq annihilation in the case of WW production). Their explicit 
expressions up to % = 3 are given in Ref. [26J in terms of the universal perturbative coefficients 
A q l \ Aq \ A q 3 \ B^' N , B^ N . In particular, the LL function #W depends on the coefficient A q , 
the NLL function g K N ' also depends on A\ and B\ and the NNLL function g K N ' also depends 
on A q 3 ^ and B^ N . All these coefficients are known except A q 3 \ In our quantitative study (see 

Sect . [3]) we assume that the value of A q ^ is the same as the one [301 EI] that appears in resummed 
calculations of soft-gluon contributions near partonic threshold. 

We now turn to the hard coefficients Ttq^ a ^pf- The first order contributions W^^jy in 
Eq. §5§ are known. The flavour off-diagonal part is process independent and in the MS scheme it 
reads: 

o/l) = <iy(l) _ I ,^(1) ] n 9l (R\ 

K q q^ gq ,N ^gg^qg,N 2 (N + 1)(N + 2) 7q9 ' N /i| ' 1 ' 

where 7^ N are the LO anomalous dimensions. The flavour diagonal coefficient is instead process 
dependent, and, as shown in Refs. [23j [21], it can be expressed in terms of the finite part A ww 
of the one- loop correction to the Born subprocess, computed in Ref. p]. We have 



n ww(i) _ c ( 1 , f \ , l A ww _ ( o(i) , ln Ml\ in Ml + 2 ~(i) l n 9l (9) 

The second order coefficients \ in Eq. ([5]) have not yet been computed. 

We finally consider the finite component in Eq. (Q. This contribution has to be evaluated 
starting from the usual perturbative truncation of the partonic cross section and subtracting the 
expansion of the resummed part at the same perturbative order. This procedure allows us to 
combine the resummed and the finite component of the partonic cross section to achieve uniform 
theoretical accuracy over the entire range of transverse momenta. Note that, since for WW 
production at pYp W 7^ only the LO result is known (VTVT-l-l parton), we can perform the 
matching at LO only. 



*We adopt here for A ww the definition of Eq. (38) of Ref. [H]. 
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In summary, the inclusion of the functions and of the coefficient in the 

resummed component, together with the evaluation of the finite component at LO, allows us to 
perform the resummation at NLL+LO accuracy. The inclusion of the function g N (still performing 
the matching at LO) allows us to reach (almost) NNLL+LO accuracy[§. The reason why we cannot 
claim full NNLL accuracy is that the NNLL contribution a^g^ in the exponent in Eq. ([6]) is of 
the same order of the terms coming from the combined effect of a\7i^ W<y2 ' > and Lg^\ which are 
not under control. For this reason, in the following Section we will mainly rely on our NLL+LO 
prediction, regarding the NNLL+LO result as an indication of the size of NNLL effects. 



3 Results 



In this Section we present numerical results for WW production in pp collisions at LHC energies. 
We compare our resummed perturbative predictions at NLL+LO and NNLL+LO accuracy with 
the NLO ones, obtained with the general purpose program MCFM j9], and with results obtained 
with the MC@NLO event generator [28] . 

To compute the WW cross section we use MRST2002 NLO densities [32] and «s evaluated 
at two- loop order. As discussed in Sect. [2j our resummed predictions depend on renormalization, 
factorization and resummation scales. Unless stated otherwise, the resummation scale is set equal 
to the invariant mass Mww of the WW pair, whereas renormalization and factorization scales 
are set to 2Mw- The latter choice allows us to exploit our unitarity constraint and to exactly 
recover the total NLO cross section when no cuts are applied. At NLO we consistently use 
(J>f = (J>r = 2Miy as default choice, whereas in MC@NLO //p and /i# are set to the default choice, 
the average transverse mass of the Ws. 

The predictions of resummation are implemented in a partonic Monte Carlo program which 
generate the full 5-body final state (lulu + 1 parton). Nonetheless, since the resummed cross 
section in Eq. ([I]) is inclusive over rapidity, we are not able to apply the usual rapidity cuts on 
the leptonfl To the purpose of the present work, we do not expect this limitation to be essential. 
Note also that, since the resummation formalism we use is valid for the inclusive production of a 
non-strongly interacting final state (WW — > lulu in the present case), we are not allowed to apply 
cuts on the accompanying jets. 

We start by considering the inclusive cross sections. Our NLL+LO (NNLL+LO) result is 115.6 
(115.5) pb, and agrees with the NLO one (116.0 pb) to better than 1%. The cross section from 
MC@NLO is instead lower, about 114.7 pb. The above difference is due to the different choice of 
the scales, and to the different convention in the choice of the electroweak couplings adopted in 
MC@NLO. 

In Fig. [1] (left) we show the corresponding pt distribution, computed at NLL+LO (dotted), 

§The reader should not be confused by this notation: the NLL+LO and NNLL+LO results include the complete 
0(as) real contribution through the finite component in Eq. ^ as well as the full virtual 0(as) correction through 
the nJgZq^N coefficient in Eq. ©. As a consequence, our NLL+LO and NNLL+LO results contain the complete 
NLO correction plus resummation of the large logarithmic contributions at p^ w <c Mww ■ 

^Note that this is not a limitation of principle: the resummation formalism can be extended to the double 
differential transverse momentum and rapidity distribution. 
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Figure 1: Left: comparison of the transverse momentum spectra of the WW pair obtained at 
NLL+LO, NNLL+LO and with MC@NLO. No cuts are applied. Right: NNLL+LO and MC@NLO 
results normalized to NLL+LO. 

NNLL+LO (solid) and with MC@NLO (dashed). The NLO result, not shown in the plot, diverges 
to +00 as p^ w — > 0. We see that the three histograms are very close to each other and show a 
peak around p^ w ~ 5 GeV. The agreement is confirmed in Fig. [I] (right) where the ratio of the 
NNLL+LO (solid) and MCONLO (dashed) results to the NLL+LO are displayed. We see that, 
apart from statistical fluctuations, the MC@NLO and NLL+LO agree within 2-3 % on average. 
On the contrary, the NNLL contribution tends to make the distribution harder and its effect 
increases at higher p?, always being below 10%. 

In order to study the perturbative uncertainties affecting our NLL+LO calculation, we have 
varied /if and hr by a factor 2 around the central value. We find that the effect of scale variations 
is rather small, of the order of ±1%, and comparable with the estimated accuracy of our numerical 
code. Similar results are obtained by varying fi F and in MC@NLO. 

The dependence of our NLL+LO results on the resummation scale is instead stronger. In 
Fig. [2] we show the NLL+LO prediction for different values of the resummation scale Q. We see 
that varying the resummation scale the effect on the pj> w spectrum is well visible and amounts to 
about ±10% at the peak. For lower (higher) values of Q the effect of the resummation is confined 
to smaller (larger) values of p^ w . Comparing Figs. [1] and [2] we see that the order of magnitude 
of the NNLL effect, partially included in our NNLL+LO prediction, is smaller than the spread in 
the NLL+LO result from resummation-scale variations. 

As in the case of Higgs production [25], we find that the choice Q = 2M WW gives (slightly) 
negative cross sections at very large p^ w . In order to define a range of variation of Q, we would 
like to avoid values that give a bad behaviour at large p^ w . For this reason in the following we 
will consider resummation-scale variations in the range Mww/^ < Q < M ww . 

We now consider the px spectra of the leptons. For each event, we classify the transverse 
momenta of the two charged leptons into their minimum and maximum values, p^mm an d Prmax- 
In Fig. Owe plot the corresponding p T spectra, computed at NLL+LO (solid), NLO (dotted) and 
with MC@NLO (dashes). We see that all the three predictions are in good agreement. Small 
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Figure 2: NLL+LO spectra for different values of the resummation scale Q. 

differences are visible in the peak of the pt distribution of the lepton with larger p^: the peak 
predicted by MC@NLO and by the NLL+LO calculation is slightly lower than the one from the 
NLO calculation. The renormalization- and factorization-scale dependence of the results, defined 
as above, is still of about ±1%. The scale uncertainty of the NLL+LO result is dominated by 
resummation-scale variations and it is about ±2 — 3% on average. 

We then examine the impact of resummation when cuts on the final state are applied. We 
start by defining the following selection criteria, taken from the study of Ref. |29j . 

Cuts A: 

• The pt of the charged leptons should be larger than 20 GeV. 

• The invariant mass mu of the charged leptons should be smaller than 80 GeV. 

• The missing pt of the event should be larger than 20 GeV. 

• The azimuthal separation A0 of the charged leptons in the transverse plane should be smaller 
than 135°. 

These cuts basically select a pair of If' s, suppressing events with lepton pairs originating from the 
decay of a Z. The corresponding NLL+LO cross section is 21.03 pb, about 5% smaller than the 
NLO result of 22.10 pb. The MC@NLO cross section is instead 21.16 pb. The differences between 
these results are not unexpected: both the NLL+LO and the MC@NLO calculations enforce a 
unitarity constraint which ensure the correct NLO normalization is recovered when total cross 
sections are considered. When, as in this case, cuts are applied, higher order effects are present 
that make NLL+LO and MC@NLO results generally different from NLO. 
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Figure 3: Distributions in p l Tmin and p l Tmax . No cuts are applied. 

In Fig. H]we show the distributions in p l Tmin and p^ max : as in Fig. [31 no significant differences 
are found between the three histograms: these distributions are still reliably predicted by the 
fixed-order NLO calculation. This conclusion is confirmed by studying scale variations. The 
effect of renormalization- and factorization-scale variations is still very small, and the impact of 
resummation-scale variations on the NLL+LO result is again of about ±2 — 3% on average. 

In the search for the Higgs boson in the H — > WW — > lulu channel an important difference 
between the signal and the background is found in the A0 distribution. Since the Higgs is a 
scalar, the charged leptons tend to be produced quite close in angle. As a consequence, the signal 
is expected to be peaked at small values of A<f>, whereas the A0 distribution for the background 
is expected to be reasonably flat. It is thus important to study the effect of resummation on this 
distribution, which is also known to be particularly sensitive to spin correlations [33] . 

In Fig. [5] the A0 distribution at NLL+LO (solid) is compared to the one at NLO (dotted) and 
from MC@NLO (dashes). The upper panel shows NLL+LO and MC@NLO results normalized to 
the NLO prediction. Note that both the NLO and the NLL+LO calculations fully include spin 
correlations, whereas MC@NLO neglects spin correlations in the finite (non factorized) part of the 
one-loop contribution. We see that the agreement between the three results is excellent, showing 
that the approximate treatment of spin correlations in MC@NLO is accurate enough. Comparing 
NLL+LO and MC@NLO predictions with NLO ones we also see that the effects of resummation 
appear negligible in this situation. As in the case of Fig. HI this conclusion is confirmed by studying 
scale variations. The effect of renormalization- and factorization-scale variations is still of about 
±1%. The effect of resummation-scale variations on the NLL+LO prediction, not shown in the 
plot, is about ±2%. 

We now consider the application of stronger selection criteria [29J , designed for the search of a 
Higgs boson with mass Mh = 165 GeV. 
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Figure 4: Distributions in p l Tmin and p^ max . Cuts A are applied. 

Cuts B 

• For each event, p l Tmin should be larger than 25 GeV and p l Tma _ x should be between 35 and 50 
GeV. 

• The invariant mass ma of the charged leptons should be smaller than 35 GeV. 

• The missing pt of the event should be larger than 20 GeV. 

• The azimuthal separation A0 of the charged leptons in the transverse plane should be smaller 
than 45°. 

• A jet veto is mimicked by imposing that the transverse momentum of the WW pair should 
be smaller than 30 GeV. This cut is perfectly legitimate in our resummed calculation and is 
exactly equivalent to a jet veto at NLO. 

These cuts further select the small A0 region. The jet veto is usually applied in order to reduce 
the it contribution, which is expected to produce large-p^ 6-jets from the decay of the top quark. 

The new cuts reduce the number of WW events by an additional factor of 35 with respect 
to cuts A. The NLL+LO (MC@NLO) accepted cross section is 0.599 pb (0.570 pb) which should 
be contrasted with the NLO result, which is 0.691 pb, about 20% higher. This relative large 
difference is due to the fact that the new cuts enhance the relevance of the small-p^^ 7 region, 
where the NLO calculation is not any more reliable. 

In Fig.[6]the p l Tinin and p l Tma _ x distributions are presented. The upper plot shows the two spectra 
for central values of the scales. We see that although the three predictions are still in reasonable 
agreement in shape, differences are now evident. In particular, the p^mm distribution at NLO is 
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Figure 5: The A<p distribution when cuts A are applied. The upper part of the plot shows the 
NLL+LO and MC@NLO results normalized to the NLO one. 

definitely steeper than the other two. Comparing NLL+LO and MC@NLO spectra, we see that 
the former are steeper than the latter: with the application of stronger cuts the differences between 
NLL+LO and MC@NLO predictions are clearly amplified. The lower plot shows the uncertainty 
band at NLO and NLL+LO. At NLO the band is obtained by varying p, F = p, R between M w and 
4M\y: the effect of scale variations is now more visible and is of the order of a few percent. The 
corresponding scale variations have marginal effect on the NLL+LO and MC@NLO results: this 
means that when the small p^ w region is selected, scale uncertainties increase at fixed order, while 
remaining small when resummation is included. At NLL+LO the band is obtained as before by 
varying the resummation scale Q between Mww and 4Mww, and the effect is of about ±3 — 4%. 

In Fig. [7] (left) the A0 distribution is displayed. The NLO and NLL+LO bands (right), 
computed as in Fig. El are also shown. We see that the shapes of the three results are still in good 
agreement with each other, although a slightly different slope of the NLL+LO result with respect 
to MC@NLO and NLO ones starts to appear. 

We finally consider the transverse-mass distribution of the WW system. We choose to define 
the transverse mass according to Ref. [34J. We start from the transverse energy of the charged 
leptons and of the neutrinos, which are expressed as 

Etu = ^P 2 tu + < $T=\ltf T + rnl v , (10) 

where p TU and ji T are the total transverse momentum of the charged leptons and the missing 
transverse momentum, respectively. If the Higgs boson mass is close to the WW threshold, the 
invariant mass of the neutrinos in the expression of L^t can be approximated by mu and we can 
define 

M TWW = y/(f T + E TU f - (pru + j* T ) 2 . (11) 
In Fig. [S] (left) we show the transverse-mass distribution defined above, computed at NLO 
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(dotted), NLL+LO (solid) and with MCONLO (dashes). The inset plot shows the NLO and 
MC@NLO results normalized to NLL+LO. The plot on the right shows the uncertainty bands 
computed as in Figs. I6TIT1 All the three histograms are peaked at M TWW between 140 and 150 
GeV. Comparing the shapes of the histograms, we see that at NLO the shape is fairly different 
with respect to NLL+LO and MC@NLO. The NLL+LO and MC@NLO distributions also show 
small differences: the NLL+LO result is steeper and softer than the MC@NLO one. 



4 Summary 

In this paper we studied the effects of soft-gluon resummation in WW pair production in hadron 
collisions. We performed a calculation that, using the helicity amplitudes computed in Ref. [7j, 
achieves uniform NLO accuracy over the whole phase space but consistently includes the resum- 
mation of the large logarithmic contributions at p T /w <C M ww . We presented predictions for the 
p T 7W spectrum at NLL+LO and (almost) NNLL+LO accuracy, and compared our results with 
those obtained with the MC@NLO event generator, finding good agreement. 

We then examined a few charged-lepton distributions: the px spectrum of the lepton with 
minimum and maximum px, the azimuthal separation of the charged leptons in the transverse 
plane and the transverse-mass distribution, when different sets of cuts are applied. The px spectra 
of the charged leptons are generally well described by the NLO calculation but the effect of 
resummation becomes visible when hard cuts are applied. The A0 distribution is important to 
discriminate a Higgs-boson signal over the background in the H — ► WW — ► lulu channel, and is 
also known to be particularly sensitive to the effect of spin correlations in the W's decay. The 
latter are fully taken into account at NLO and in our NLL+LO calculation, whereas MC@NLO 
neglects spin correlations in the finite (non-factorized) part of the one-loop contribution. All the 
predictions for the A0 distribution agree well and resummation effects appear small even when 
hard cuts are applied. We also examined the transverse mass distribution of the WW pair. The 
shape of the NLO distribution differs substantially from NLL+LO and MC@NLO. On the other 
hand, the NLL+LO and MC@NLO results are still in reasonable agreement. 

We finally note that the resummed predictions presented in this paper were obtained in a 
purely perturbative framework. Two kinds of non-perturbative effects may be considered: those 
due to hadronization, and those due to an intrinsic transverse momentum (intrinsic px) of the 
partons in the colliding hadrons. 

In the present paper we were mainly concerned with leptonic observables. As a consequence, 
we do not expect hadronization effects to be particularly relevant. This expectation is supported 
by the good agreement we find between resummed results and MC@NLO predictions, in which 
hadronization effects are taken into account. 

As far as intrinsic-pr effects are concerned, it is known (see e.g. Ref. [16] and references 
therein) that transverse-momentum distributions are affected by such kind of non-perturbative 
effects, particularly at small transverse momenta. These effects, associated to the large-6 region in 
impact parameter, are not taken into account in our calculation. As noted in Ref. [25] , these non- 
perturbative effects have the same qualitative impact as the inclusion of higher-order logarithmic 
contributions, i.e., they tend to make the pt distribution harder. The precise quantification of 
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non-perturbative effects in WW production is left for future investigations. 
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Figure 6: Distributions in p Tmin and p TmSLX when cuts B are applied. In the lower plot the pertur- 
bative uncertainty bands at NLO and NLL+LO are shown. 
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Figure 7: Left: A(j) distribution when cuts B are applied. Right: the corresponding NLL+LO and 
NLO uncertainty bands are shown. 
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Figure 8: Left: Mtww distribution when cuts B are applied. The inset plot shows NLO and 
MC@NLO results normalized to NLL+LO. Right: NLL+LO and NLO uncertainty bands computed 
as in Fig. [7| 
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